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Introduction 

Although there seems to be no commonly accepted definition of integrability for general 
systems of partial differential equations (PDEs), 121113011 , some features of these systems are 
beyond any doubt. These are: 

- infinite hierarchies of (possibly, nonlocal) symmetries and/or conservation laws, 

- (bi-)Hamiltonian structures, 

- recursion operators. 

Of course, these properties are closely interrelated, even if this is not so straightforward as 
it may seem at first glance (see, e.g., |21). 

Methods to establish the integrability property are numerous and often reduce to finding 
a Lax pair for a particular PDE. In spite of their efficiency, all these methods possess two 
drawbacks that are quite serious, to our opinion: 
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- they are unable to deal with general PDEs directly and to be applied require reduction 
of the initial system to evolution form (which is not always done in a correct way, see, 
e.g., the discussion in |7|), 

- they lack a well-defined language for working with nonlocal objects that are essential in 
most applications. 

In a series of papers l7l[8l llO|[TT]| 1211 1311 1411 151 we presented our approach to integrable 
structures which is, at least partially, free from these disadvantages. In I16II a review of the 
geometrical theory illustrated with many examples from the above papers can be found. 

Essentially, the approach is based on two constructions naturally associated to the equa- 
tion at hand: those of tangent and cotangent coverings (see Sec.|5]l that serve as exact concep- 
tual counterpart to tangent and cotangent bundles in classical differential geometry. These 
coverings are also differential equations, and we prove that all operators responsible for the 
integrability properties (i.e., Hamiltonian, symplectic, and recursion ones) can be identi- 
fied with higher or generalized symmetries and cosymmetries of these equations (Sec. |4]i. 
Consequently, all computations boil down to solving two linear equations: 



where ££■ is the linearization operator of the initial equation and £*c, is its formally adjoint 
(see Sec.|2]and[3]l lifted to tangent or cotangent coverings. The solutions that possess certain 
additional properties (expressed in terms of the Schouten and Nijenhuis brackets and the 
Euler operator) deliver the needed operators. 

In this paper, we tried to expose our method in such a way that the interested reader 
could use it as an operational manual for computations. As a tutorial example we chose the 
Korteweg-de Vries equation because all results for this equation are well known and can 
be easily checked, even by hand. More complicated equations are also briefly discussed. 
In particular, we will consider equations in more than two independent variables which 
are naturally presented in non-evolutionary form (Kadomtsev-Petviashvily and Plebanski 
equations). The related examples appear here for the first time. In Sec. [TO] we describe a 
freely available software package that was used in our computations 1271 . 

1 Differential equations and solutions 

Our general framework is based on jT], but we will give a coordinate-based exposition here. 

For an equation in n independent variables x' and m unknown functions we consider 
the jet space J°°{n,m) with the coordinates x' , m^, where ; = 1, . . . ,n, _/ = 1, . . . ,m and a = 
iyii . ■ where <iii < n and ii < !2 < • ■ • < is an ordered multi-index of finite, but 
unlimited length | ct | . Denote by ;r : 7°° (n, m) — > R" the projection to the space of independent 
variables x^ , . . ., x". 

The vector fields 



and 




u 



d 



(1) 



where ai is the ordered multiindex obtained by ct and / after reordering, are called total 
derivatives. They define an «-dimensional distribution on J°°{n,m) that is called the Cartan 
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(or higher contact) distribution and is denoted by "rf. Dually, the Cartan distribution is the 
annihilator of the system of differential 1 -forms 



CO, 



\iy = du'a -Y^u'„;dx\ j=l,...,m, \a\>0, (2) 



the so-called Cartan, or higher contact forms. 
Let a differential equation be given by 



^(•■■'^''■••'^'■■•j=^' '=1'---'-- (3) 

Then we consider all its differential consequences, or prolongations, 

D„f'=0, 1=1,. ..,r, \a\>0, (4) 

where Da = Di^ o ■ ■ ■ o/),^ for o = ii . . .i/^. We take the (usually, infinite-dimensional) hy- 
persurface <£' C J°°{n,m) defined by zeros of (|4li as a geometrical image of 

We assume that the following two conditions hold: first, if G is a function on J'^{n,m) 
and G\g = then Gis a differential consequence of (F', . • • ,F'') \ second, for any differential 
operator in total derivatives A such that A(F') = 0, / = 1, . . . ,r, we have A\g = Q. The 
second condition excludes gauge equations from our consideration. 

Differential operators in total derivatives are said to be -differential operators. 

By lIU, all total derivatives are tangent to S and thus the Cartan distribution induces an 
w-dimensional distribution on S, which we call by the same name. Its «-dimensional integral 
manifolds are solutions of S. 

Example 1 Consider the KdV equation 

U, = UUjc + U^xx- (5) 

The corresponding jet space is /°°(2, 1) with the coordinates 

limes k times 

and the total derivatives 

ax .J^Q dui^k at ^.^q dui^t 

Thus, the space S is defined by the relations 

Ui,k+\ =0^-Of(Mo.OMl,0-|-M3,o), 

from where it follows that all partial derivatives containing t may be expressed via the ones 
containing x only. The elements of the first group are called principal derivatives, and the 
remaining ones are called parametric derivatives L20ll . 

The functions x, t, and the parametric derivatives m, = m, rnay be taken for coordinates 
on (internal coordinates). In terms of these coordinates, the total derivatives are 

Oa = ^ + ^m/+i^, a = 4- + L^i(""i +"3)^> (6) 
dx dui dt dUi 

while the Cartan forms on S are given by 

(Oi = dui — M,+ i dx — D'^(uu\ +U2)dt, i>0. 
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Example 2 The KP equation 



Uyy — Utx 11^ ^l^xx J 2 ^xxxx (7) 

is defined in the jet space 7°°(3, 1). Its total derivatives can be written as 

(9 (9^ d ^ d d 

o/f dua dx dua dy dua 

where P is the set of multi-indexes that do not contain more than one instance of y. This is 
the set of multi-indexes of all parametric derivatives. Note that in D,. the coordinates M(jy, 
where ct £ P is of the form a = Ty, are principal, and must be replaced with 

n 2 1 

i^T I Utx ~U-x~ UUxx — — ^xxxx 



1 



2 Linearization and symmetries 



A (non-trivial) symmetry of the distribution ^ (on J°°{n,m) or on (a) is a ;r-vertical vector 
field that preserves this distribution. In other words. 



diiL 



is a symmetry if and only if 



[Di,X] = 0, 



(8) 



Equation ^ implies that 



1,. 



|ct|>0, 



from which it follows that any symmetry is of the form 



(9) 



on J'"[n,m). Here (p = (9^, .... (p'") is an arbitrary smooth vector function onJ°°{n,m). Vec- 
tor fields of the form ^ are called evolutionary and 9 is called the generating function (or 
characteristic). Generating functions are often identified with the corresponding symmetries 
and form a Lie R-algebra with respect to the Jacobi (or higher contact) bracket 



(10) 



Given an equation S, its symmetry algebra symd" is formed by the generating func- 
tions (p, satisfying the system of linear equations (defining equations for symmetries) 



(11) 



where the operator £g is defined as follows. For any smooth vector function F = {F^ . . ,F'') 
on J°°(n,m) set 





dF" 


E 

a 


diVf, 
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This operator is called the linearization of F. It is a matrix differential operator in total 
derivatives acting on generating functions and producing vector functions of the same di- 
mension as F. Then, if ra is defined by Eq. ^ one has 

The restriction of £f to (f is well defined since £f is an operator in total derivatives. 
Example 3 For the KdV equation ([5]l, the defining equation for symmetries is 

Di{(p) = ui(p + uDA(p)+Dl{(p), (12) 
where (p = (p{x,t,u,u\, . . . ,ui^) and D^, D, are given by 

3 Conservation laws and cosymmetries 

Consider an equation ra C J°°{n,m) that, for simplicity, we will assume to be topologically 
trivial. This means that we assume that the equation ^ and its prolongations (|4]i be topolog- 
ically trivial. This holds at least locally for most equations, and it is a reasonable hypothesis 
in computations. Define the spaces A,^((?') of horizontal dijferential forms that are generated 
by the elements 

a = adx'^ A---Adx''', aeC°°{S'). 
The horizontal de Rham differential : A^^ — > A^,^^^ is defined by 

dhCO = 52 ^'(«) dx' Adx'' A---A dx''' , 
i 

and di, o di, = 0. 

An [n — l)-form (O 6 A^"^' [S) is a conservation law of S if J;,a) = 0. A conservation 
law is trivial if it is of the form dh 6 for some e eAll-^iS-). Thus we define the group of 
equivalence classes of conservation laws to be the cohomology group CL[S') = [S) = 
kerdh/ imdh- 

To facilitate computation of nontrivial conservation laws, their generating functions (or 
characteristics) are very useful |1, p. 29] (see also I24II25II26I ). To present the definition, 
recall that for any "^-differential operator A: V — > W its formal adjoint A* : W ^ V i& 
uniquely defined, where A = hom(A,Ay"). Namely, for a scalar operator A = Y.a'^cy^a one 
has 

A*=Y^{-\t\D^oaa 
a 

and for a matrix one A = \\Aij\\ we set 

A* = ||A*.||. 

Now, let V be the space to which F = (F f"'') belongs, F-' being the functions from 
equation Consider a conservation law a ^A'^^^{S) and extend it to a horizontal (n — 1)- 
form (B on J°°{n,m). Then, due to the regularity condition, since dij{a)) must vanish on tf, 
we have 

dh{a))=AiF) 

for some "^-differential operator A : V — ^ A/"(/°°). The characteristic Xj/a, of a> is defined by 

V/« = A*(1)|^. (13) 
Then a> is trivial if and only if \j/(o = 0. 
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Example 4 Let S he, & system of evolution equations 

«/ = f^[x,t,...u'j^,...), j,a= l,...,m, 

where subscript k indicates the number of derivatives over x. Let CO = X dx + T dt be a 
conservation law. Then its characteristic is 

- ^\ 

where 

are the variational derivatives. 

Generating functions satisfy the equation 

e*A¥) = o- (14) 

Arbitrary solutions of this equation are called cosymmetries of S. The space of cosymme- 
tries is denoted by cosym((f ). Thus, construction ( I13l > determines the inclusion 

5 : CL(<r) ^cosym(<r). (15) 

Example 5 For the KdV equation ([5J its cosymmetries are the solutions of 

A(v/) = mD,(v/)+£».?(V^), (16) 
as it follows from the expression of the adjoint operator of ig (see (I12t). 

4 Nonlocal extensions 

An invariant, geometric way to introduce nonlocal objects (e.g., variables, operators) to the 
initial setting is based on the concept of differential coverings 1 1 , Chapter 6], see also the 
collection of papers 11171 . To introduce the concept, consider two well known examples. 

Example 6 Let 

U, = UUjc + U^x 

be the the Burgers equation. Then the Cole-Hopf substitution 

u - - — 

V 

introduces a new variable v which is nonlocally expressed in terms of the old one and satis- 
fies the heat equation 
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Example 7 The Miura transformation 



1 2 
Wx = U + -W 



introduces a nonlocal variable w and transforms the KdV equation ^ to 

1 2 

Wt =Wxxx--zW Wx, 

o 

i.e., to the modified KdV equation. 

The geometrical construction that generalizes these and similar examples is as follows. 
Consider two equations S' C J°°{n,m) and S C J°°(n,in) with the same collection of inde- 
pendent variables jc' , . . . , jc". A smooth surjection T : ^ ^ ^ is called a dijferential covering 
(or simply a covering) of tahy S if its differential takes the Cartan distribution on S to that 
on S. Coordinates in the fiber of T are called nonlocal variables. 

In coordinates, this definition means that the total derivatives on S are of the form 

bi =Di+Xi, i= l,...,w, 
where X,- are T-vertical vector fields that satisfy the relations 

Di{Xj) - Dj{Xi) + [Xi,Xj] =0, l<i<j<n. (17) 

Any "^-differential operator A on (a can be lifted to a "^-differential operator A on £' sub- 
stituting all total derivatives D, by 6, . 

Example 8 Consider the case of three independent variables x, y, z- Let CO = adx Ady + 
bdxAdz be a conservation law of ^ which has only two nontrivial components. Set £" = 

<r X K°° and 

Dx = Dx +Wx^+ Wxx 1 , 

aw awx 

d , , d 

Dy = Dy + a— +Dx{a) - — + • • • , 
aw aWx 

D, = D, + b^ +Dx{b)^ + 
aw aWx 

where w, Wx, Wxx, ... are coordinates in Then T : x M°° — is a covering. 

The same construction works for n / 3. If n = 2 we get a 1-dimensional covering. 
For w > 2 it is important to remember that the initial conservation law must have only two 
nonzero components. 

An interesting and a well known class of coverings comes from the Wahlquist-Estabrook 
construction I28II29I . 

Example 9 For the KdV equation, consider the coverings 

Dx = Dx+X, D,=D,+T 

such that the coefficients of the fields X and T depend on mq, "i, "2 and nonlocal variables 
only. Then one can show that for any such a covering 

X = UqA + uqB + C, 
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T = {2UQU2 -u\ + ^ul^A+{u2 + ]^ui^B + ui[B,C] + ]^[B, [C,B]] + mo[C, [C,S]] +D, 

where the fields A, B, C, and D are independent of if, and fulfill the commutator relations 

[A,B] = [A,C] = [C,D] =0, 
[B,[B,[B,C]]] = Q, [5,D] + [C,[C,[C,S]]] = 0, (18) 

[A,D] + i[C,5] + ^[S,[C,[C,5]]] = 0. 

Denote by g the free Lie algebra with four generators and defining relations ( 118b . Then any 
covering satisfying the above formulated ansatz is defined by a representation of q in the Lie 
algebra of vector fields on the covering. 

The next example is crucial for the subsequent constructions. 

Example 10 The class of coverings that will be introduced in this example delivers the so- 
lution to the following factorization problem. Assume that two "^f -differential operators 

A:V^W, A':V'^W' 

are given, where V, W, V', W' are spaces of smooth vector functions on of dimensions /, 
s, /', and s', respectively. The problem is: find all operators A: y — > V' such that there exists 
a commutative diagram 

V — ^ W 

i [' 

V > W 

A' 

with some operator S: W — > W'. Since the problem has obvious trivial solutions of the 
form A=A'oA,'we are interested in the space 

{A\A'oA=BoA} 
^(4,4' =^^-^ T-^ (19) 

^ ' ' {A\A=A'oA} 

consisting of nontrivial solutions. Obviously, elements of , A'} define maps from kerzi 
to kerA'. 

Assume that ^ C J°°(n,m) with the coordinates jc', u". Consider the space J°°{n,in + 1} 
with additional coordinates v^. There is a natural surjection ^ : J°°{n,in + l) J°'{n,m), 
and let ra^ be defined by 

4(v) =0 

in ((?). Then ^: ra^ ^ (a is evidently a covering. We call it the A-covering. 
To every operator 

a 

there corresponds a vector function 



A 



a.a 
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on Sa ■ The function 4>4 vanishes if and only if A = A' o A for some operator A', so that the 
function <Pa corresponds to the equivalence class of A. 

As was indicated above, the operator A' can be lifted to A' in the covering under con- 
sideration. Consider a function <i>A that satisfies the equation 

A'(4>a)=0. 

It can be shown 1101 that such functions are in one-to-one correspondence with the elements 
of i/(A , A') given by ( fT9] >. 



5 Tangent and cotangent coverings 

These two coverings serve, in the geometry of differential equations, as counterparts to 
tangent and cotangent bundles in differential geometry of finite-dimensional manifolds II16I . 
Technically, both can be obtained as particular cases of ExamplellOlabove. 

Consider an equation S C J°°{n,m) and take the operator for A in the construction of 
the A -covering. In more detail, this means that we consider a new dependent variable q = 
{q^,...,q"^) and extend S" by the equation 

The resulting equation is denoted by S^{S) C J"'{n,2m) and the projection t: ^^{S) — > S 
is called the tangent covering to S. A important property of T is that its holonomic sections, 
i.e., sections (p: S ^ 3^{S) that preserve the Cartan distributions, are symmetries of S. 

Example 11 The tangent covering for the KdV equation is the system 

qt = Uxq + uqx + qxxx (20) 

that projects to S hy {x,t,...,Uk,...,qi,...) {x,t,...,Uk,...). 

Dually, consider the jet space J°°{n,m + r), where r is from with the new dependent 
variable p = (p' , . . . , p'') and extend by 

e*Ap) = o- 

The equation ) C J°°{n,in + r) together with the natural projection T* : ^*(tf ) S 

is the cotangent covering to S. Its holonomic sections are cosymmetries of S. Note also 

that ^*{S') is always an Euler-Lagrange equation with the Lagrangian L = F^p^ -\ h 

F'p''. 

Example 12 If S is the KdV equation then <T* [S) is of the form 

Ut = UUx + Uxxx J 

Pt = UPx+Pxxx, (21) 



cf. (Ull. 
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There exist two canonical inclusions, 

v: sym(<f) ^CL(5^*(r) 

and 

V*: cosym(^) ^ CL(.5^<r) 

(defined below), that play a very important role in practical computations described in 
Sec.[6H9]below. 

Moreover, one can prove (see 1161 ) that the image of v consists of the elements that 
are fiber-wise linear with respect to the projection t: — > S". Dually, the image of v* 
coincides with conservation laws that are linear along the fibers of the tangent covering 

In what follows, we assume that r (the number of equations) equals m (the number of 
unknown functions). 

Consider a symmetry (p of the equation (f and extend it to some (p defined on the ambient 
space J°°{n,m + r). Then due to the Green formula jT] p. 191] (see 1241 p. 41] for more 
details) one has 

where (■, ■) is the natural pairing between the spaces A and A, while the correspondence (p i— )• 
03^ is a "^-differential operator. We set 

v{(p) = ®^|^.^, 

and this is a well defined operation. Note that in the right-hand side of the above formula 
and in similar situations, to simplify notation, when writing a differential form we actually 
mean its de Rham cohomology class. 

Example 13 Let (p be a symmetry of the KdV equation. Then, due to ( 112b and illi . one has 

{D,{(p) - Dl{(p) - iiD^{(p) - u^(p)p - (p{upjc + Pxxx - Pt) 

= D,{(pp)-D,{up(p - 3p,A{(p)+Dl{p(p)), 

i.e., 

v{(p) = (ppdx+ {iip(p-3p,,D,{(p)+Dl{p(p))dt. (22) 

In a dual way, if y/ is a cosymmetry of S" then extend it to some \j/ on J°°{n, 2m), use the 
Green formula 

{tp{f),q) - {\j/jF{q)) = dhCOf, 

and set 

with similar considerations as v. 

Example 14 For a cosymmetry i// of the KdV equation, using equations i ll6b and i l20b . we 
arrive at the correspondence 

= \j/qdx+ [ucjxj/ -3qxD_^{\j/) +D^.{q\l/)) dt, 

which in this particular case coincides with ll22t . 

' To be more precise, the images of V and u* consist of equivalence classes that contain at least one 
fiber-wise linear element. 
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Example 15 Consider the Plebanski (or second heavenly) equation 

Uxz + M/y + UnUxx - "u- = 0. 



(23) 



Its linearization is self-adjoint (so that it is Lagrangian), hence the tangent covering coincides 
with the cotangent one; they are defined through the equation 

Pxz. + Ply + UtiPxx + UxxPti - 2utxP,x = 0. 

The M- translation symmetry d/du has the generating function (p = 1, so we have 

Pxz + Ply + UltPxx + UxxPtt - ^UtxPix 

= DxiPz + UttPx - U,xPt)+D,{py + UxxPt - UlxPx) = 0, 

thus the corresponding conservation law is 

1^ ( 1 ) = (Py + I'xxPi — UixPx) dxAdyAdz+ {u,xPt — Pz — u„ Px) dt Ady Adz. (24) 



6 Recursion operators for symmetries 

Consider an equation S" and its tangent covering T : — >■ ^. Let be the lifting of the 
linearization operator to this covering and 4" = (4>' , . . . , 4"'"), where 



is a solution of the equation 



Then, by Example [TOl the operator 



a' a 



4(4>)=0. 



takes elements of keri^ to themselves. In other words, these operators act on symmetries 
of (a and thus are said to be recursion operators. 

Example 16 Let 

Ut = Uxx 

be the heat equation. Then the space of the tangent covering is given by 

Ut = Uxx, Ql — Qxx 

and the defining equation for recursion operators is 

b,{^)=Dl{<P), (25) 

where 

,>0 ,>0 
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(recall that qi = Qx-.-x, i times x) and 

<t> = cPq + a^ q\-\ — ■+a^qi(, a;6C~(^). 
The first three nontrivial solutions of i l25l l are 

X 

<i>m = q, 'i'm = q\, 'i>n=tqi + -, 

to which there correspond the recursion operators 

Mm = id Mw= D„ i = fD., + ^ . 

It can be shown that the entire algebra of recursion operators for the heat equation is gener- 
ated by Si\Q and M\\ with the only relation 

[.^io,.^ii] = i 

This algebra is isomorphic to the universal enveloping algebra for the 3-dimensional Heisen- 
berg algebra. 

Example 17 Taking the tangent covering | |2QI ) to the KdV equation and solving lg{^) = 
in this covering leads to the only solution <P = q, i.e., it yields the identity operator. This is 
not a surprise, because it is well known that the KdV equation does not admit local recursion 
operators. 

Let us now take a cosymmetry \j/ = I of the KdV equation and consider the correspond- 
ing conservation law on the tangent covering 

V*{1) = qdx+ {uq + q2)dt. 

Then the corresponding nonlocal variable (see Example[8ll satisfies the equations 

In the extended setting the linearized equation (.g{<i>) = 0, where 

<i>=AiQ^ +a°9 + a''?l H ^a'^qt, 

acquires a new nontrivial solution 

2 1 1 
=^'2 + 2"?+ 3"iG 

to which the Lenard recursion operator 

,21 , 
M = Dl + -u+-mD^^ (26) 

corresponds. 

If we consider the next cosymmetry \j/ = u then the corresponding nonlocal variable will 
be defined by 

G.v = qu, = qu^ + q^^u - q^u^ + qUxx, 

which will lead to another recursion operator 

D'l + ^uDI+2uiD, + ^{u^ + 3u2) + ^{uui +U3)D;'^ + ^uiD'^ou, 
the square of ( 126b . 
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Given an equation S and two recursion operators .^i, .'%2 ■ sym(o sym^, the Nijen- 
huis bracket 

[[^1 ,^2]] : symS" x symS" symS' 

is defined by 

[[^l,^2]](<Pl,<P2) = {^l(<Pl),i^2(<P2)} + {^2(<Pl),^l(<i02)} 

-^l({^2(<Pl),<P2} + {<Pl,^2(<P2)})-^2({^l(<Pl),<P2} + {<Pl,-^l('P2)}) 

+ (,'^10^2+^20^1){(P1,(P2}, 

where {•, •} is the Jacobi bracket given by Equation dlOl l. An operator M is called hereditary 
if [[^,^]] = 0. This property was introduced in (3]|4]|5l|6l. It can be shown that if this prop- 
erty holds and M is invariant with respect to a symmetry (p then the symmetries (p, =M'{(p) 
form a commuting hierarchy (see, for example, |,9i Ch. 4, Sec. 3.1]). 



7 Symplectic structures 

Let us now solve the equation 

on the tangent covering for f of the form 



The corresponding '^-differential operator 



due to the results of Example [TO] takes symmetries of S" to its cosymmetries. 
Let us call a conservation law CO admissible (with respect to y) if 

5{co) e imy, 

where the operator 5 : Ch{^) — > cosym((o ) is from l ll5t . Then ^ determines a bracket on 
the set of admissible conservation laws by 

{co,a'}.y = La^{c'), (27) 

where Lg^{c') = £c'{<p), <p is such that 5{co) = and c' is a representative of the coho- 

mology class O)'. This bracket is well defined on equivalence classes of nontrivial admissible 
conservation laws. It is skew-symmetric if 

5^*oig=tgoy, (28) 



i.e., if f J is a self-adjoint operator. In the case of evolution equations, ( I28t means the 5^ 
is a skew-adjoint operator. 

Let now an operator satisfy Eq. I l28t . Let us deduce the conditions under which 
the bracket I l27t fulfills the Jacobi identity. To this end, consider a vector function 9 = 
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(<p', . . . , 9'") on/°°(n,m). Then, since S satisfies the regularity condition from Sec.[T] Equal- 
ity ( 128b implies 

where A,f, is a "^-differential operator. Put A(p = A(p\^ and define the map 

5J^: sym(o X symi^ — >■ cosymd" 

by 

(5^)(<Pi,<P2) = {dn'^){(P2) - {d<^y){(pi) + A*^{(pi). 
Then the equality 

5^ = (29) 

guarantees the Jacobi identity for the bracket {•, - jy. Operators satisfying | |28| | and ( 129b are 
called symplectic structures. 

Example 18 The KdV equation does not possess local symplectic structures, but if we ex- 
tend the space S^S with nonlocal variables and (see Section[6]l and solve the equation 

4(f) = 

for f of the form 

=A3Q^ +AxQ.^ +a^q + a^qi+---+a''qk, 
there will be two nontrivial solutions 



with the corresponding nonlocal symplectic operators 
and 



8 Hamiltonian structures 



Consider now the cotangent covering to S and the lifting j 
this covering. If ^ = {<i>^ . . , <P"') is of the form 



a'a 



of the linearization operator to 



and satisfies 
then, due to Example [lOl the operator 



takes cosymmetries of the equation (a to its symmetries. 
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For any such an operator, consider the bracket 

{(B,(B'}.^=L^.(5«)(ffl') (30) 

defined on the space CL{^) of conservation laws. Here 5 : CL{<a) — > cosym((l') is the 
map ( 115b . This bracket is skew-symmetric if 

{igoM'y = egojt. (31) 

For evolution equations one usually requires an additional condition of existence of a 
conservation law X such that m, = Jif{5x)- For such equations Oll l implies Jf* = —Jif. 

For any two "^-differential operators Jffi, M^'- cosym(<^') — > sym(<l') satisfying ( 13 lb 
define their Schouten bracket 1161 p. 13] (see 1241 p. 226] for a more general definition) 

: cosym(fr) x cosym((g') sym((?) 

by 

[[J^,^2]](V^1, V^2) = ^(L^^v'.(V'2)) -^2(-L^^,V., (V^2)) 

+ {,i«1(V/2),^2(V^l)}-{^(Vl),^2(V2)}, (32) 

where {■, •} is the Jacobi bracket (TJ p. 48] and 

for any 9 6 sym((o') and \\f 6 cosym((o'). 

An operator ,3^ satisfying ( 131b and such that 

=0 

is called a Hamiltonian structure on S . For such a structure, the bracket ( 130b satisfies the 
Jacobi identity. Two Hamiltonian structures J^, ^^2, are compatible if 

[[^,^2]]=0. 

Equations admitting compatible structures are called bi-Hamiltonian. One can apply the 
Magri scheme I19II to generate infinite series of commuting symmetries and conservation 
laws. 

Example 79 In the case of the KdV equation, we solve the equation 

5(4)) = mi4>-FmA,(4>)-f£>^(4>), = aV + aVi + 

in the cotangent covering (Example 1 12b and obtain two nontrivial solutions 

^21 

9\=p\., 93= P3 + :^upi + -^uipQ 

with the corresponding (and well known!) Hamiltonian operators 
Ml=D„ M,=dI+'^uD,+ ^ux. 
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Consider the ^-translation which is a symmetry of the KdV equation and the corre- 
sponding nonlocal variable P' in the cotangent covering. We have 

PI = pUx, = p{uUl +Ui)+p2Ul -p\U2. 

One obtains a new solution 

4 4 4 1,1 

*'5 = P5 + i^upT, +2UXP2 + -{u +'iU2)p\ + [-uui + -uj,)p- -^uxPx 

in the extended setting to which the nonlocal Hamiltonian structure 

Ms=dI + juDI+2uiD^ + ^{u^ +3u2)D;, + {^uui + ^ms) - ^uiD^^ oui 
corresponds. 

Example 20 Consider the KP equation (Example |2]i. The cotangent covering is defined 
through the further equation 

1 

Pyy — Ptx ^Pxx Y2 P^^^^ ' 

We will solve the equation 

4(4>) = D,,{4>) -D„(4>) + 2u,b,{<P) + M,,4> + i(D,,(4>) + -j^D,,,,(4>) = 0, 

with <i> = Y,a o.'^ Pa (here the sum extends to a finite number of multi-indexes of parametric 
coordinates a E P, and a'^ are smooth functions on (f ). The nontrivial solution 

4" = Pxx 

can be found. A simple computation shows that the corresponding operator = satis- 
fies Olb and fulfills the condition [[J^, ,^]] = 0, hence it is Hamiltonian. 
hi 1181 the following evolution form of the KP equation is introduced: 

Yx = My, Mx = V, Vx = W, W,, = 12(m, - UV - Xjfy) . (33) 

Then, one local Hamiltonian operator (Eq. (51) in 1181 ) is obtained. After applying the 
change of coordinate formulae for Hamiltonian operators (see I13II ) it is readily proved that 
it corresponds to the above (p. 



9 Recursion operators for cosymmetries 

Finally, let us consider the equation ^J(f ) = lifted to the cotangent covering. An oper- 
ator f that fulfills the equation takes an element \j/ £ cosym(tf ) into an element ^Pi^j/) £ 
cosym((o'), i.e., it is a recursion operators for cosymmetries. 
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Example 21 If S is the KdV equation then the equation 

possesses no nontrivial solution. But when we extend 3"* S by the nonlocal variable P' (see 
the previous section) a solution 

2 1 , 

arises to which the recursion operator 

M = Dl + -u-D^^ oui 

corresponds. 

Example 22 Consider the Plebanski equation (Example 1 15b. The linearization of this equa- 
tion is self-adjoint, hence the operators that we defined in this paper (recursion for symme- 
tries, symplectic, Hamiltonian, recursion for cosymmetries) act in the same spaces. Now, 
extend the (co)tangent covering with the nonlocal variable R defined by 

Rl = U,xPt -Pz- Utt Px, Rx = Py + UxxPt - UtxPx- 

The extension of the equation ig{<P) = Q admits the solution linear with respect to pa and 

^=R, (34) 

besides the trivial solution <Pi = p. 

In II22I the following evolution form of Plebanski equation is introduced: 

1 2 

u,=q, q, = {q^ - qy - M^.) . (35) 

^xx 

Then, a local Hamiltonian operator (Eq. (11) in [22 ]) and a nonlocal Hamiltonian operator 
(Eq. (27) in |22 |) are obtained. After applying the change of coordinate formulae for Hamil- 
tonian operators (see 1 13 1) it is readily proved that they correspond to 4>i and 4>2. Different 
changes of coordinates from ( 123b to ( 135b transform 4>i and 4>2 to recursion operators for 
symmetries or cosymmetries and symplectic operators. 

It is interesting to remark that, while in the evolutionary form all such operators are 
mutually different because they act in different spaces, in the original formulation of the 
equation they coincide with 4>i and 4>2, whose expressions are considerably simpler than 
that of their 'evolutionary' counterparts. 

10 Computer support 

The computations in this paper were done by the help of CDIFF, a symbolic computation 
package devoted to computations in the geometry of DEs and developed by P.K.H. Gragert, 
P.H.M. Kersten, G.F. Post and G.H.M. Roelofs at the University of Twente, The Nether- 
lands, with latest additions by one of us (R.V.). CDIFF runs in the computer algebra system 
REDUCE, which recently has become free software and can be downloaded here | 23 1. 

The 'Twente' part of CDIFF package is included in the official REDUCE distribution, 
and most of it is documented. Additional software, especially geared towards computations 
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for differential equations with more than two independent variables, can be downloaded at 
the Geometry of Differential Equations website http : //gdeq. org The software includes a 
user guide |27] and many example programs which cover most examples presented in this 
paper. 

Acknowledgments. We thank the referees for useful comments that helped us to im- 
prove the text. 
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